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ABSTRACT 

The  problem  of  wave  dispersion  and  stability  for  a  class  of  hot, 
inhomogeneous,  collisionless  magnetoplasmas  is  reduced  to  the  solution 
of  an  integral  equation  with  well-behaved  kernel.  Admissible  config¬ 
urations  include  those  for  which  the  externally  applied  and  internal 
ambipolar  fields  form  a  generalized  harmonic  oscillator.  The  full  set 
of  Maxwell's  equations  is  used  to  arrive  at  self-consistent  pertur¬ 
bation  fields  in  terms  of  the  equilibrium  particle  distributions.  An 
illustrative  example  treats  a  magnetoplasma  column  with  Gaussian 
radial  profile  and  Maxwellian  velocity  distribution  in  a  state  of  quasi¬ 
equilibrium. 
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I .  INTRODUCTION 


The  propagation  and  evolution  of  small  disturbances  in  a  hot, 
magnetized,  inhomogeneous  plasma  is  a  matter  of  considerable  theoret¬ 
ical  and  practical  interest  with  regard  to  both  the  wave  supporting 
capabilities  of  the  medium  and  its  stability  against  small  perturbations. 
Although  the  procedure  leading  to  descriptions  and  predictions  of  the 
response  of  the  plasma  to  small  signals  can  be  formulated  quite  explic¬ 
itly,  the  problem  has  long  resisted  definitive  analysis,  by  virtue  of 
what  appears  to  be  an  intrinsic  incompatibility  with  the  ideal  conditions 
under  which  standard  analyses  and  perturbation  techniques  are  effective. 
The  major  complicating  features  are  the  nonlinearity  of  the  system 
dynamics  and  the  inhomogeneity  of  the  medium,  which  makes  the  system 
sensitive,  via  particle  transport,  to  global  as  well  as  local  conditions. 

The  general  analysis  procedure  can  be  decomposed  into  a  number  of 
distinct  subproblems,  each  of  which  isolates  a  particular  feature  of 
the  system,  such  as  inhomogeneity,  nonlinearity,  anisotropy,  statistical 
distribution,  electromagnetic  interaction.  These  may  be  dealt  with 
individually  at  first,  then  combined  in  a  self-consistent  manner.  The 
concatenation  of  the  subproblems  links  the  tractability  of  one  quite 
strongly  to  that  of  the  preceding  one,  with  the  result  that  the 
complexity  of  the  mathematical  description  forces  recourse  to  numerical 
analysis  at  a  very  early  stage  of  the  calculation.  It  soon  becomes 
hopeless  to  direct  and  interpret  a  realistic  computation  meaningfully. 

It  is  the  purpose  of  this  paper  to  present  an  approach  to  the  problem 
which  can  be  pushed  analytically  nearly  to  completion,  leaving  to 
digital  computers  only  some  relatively  straightforward  quadratures. 

The  system  dealt  with  is,  in  the  general  formulation,  a  collection 
of  charged  particles  of  two  species,  inhomogeneous ly  distributed  in 
space,  with  arbitrary  velocity  distributions,  immersed  in  a  uniform 
and  constant  applied  magnetic  field,  in  the  absence  of  collisions. 

The  evolution  of  arbitrary  weak  disturbances  of  an  equilibrium  state 
is  sought,  for  the  determination  of  stability  and  wave  dispersion.  The 
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problem  ir  three-dimensional  and  is  treated  nonrelativistically ,  hut 
with  the  full  set  of  Maxwell's  equations  and  retaining  the  r.f.  magnetic 
field  interaction.  A  particular  class  of  plasma  configurations  is 

considered,  whose  realizability  is  discussed  and  justified  on  physical 
grounds . 

Experiments  on  inhomogeneous  magnetoplasmas  were  reported  and 
interpreted  by  Buchsbaum  and  Hasegawaf 1 1 2)  They  obtained  a  fourth- 
order  differential  equation  for  the  potential,  without  accounting  for 
the  ami, i polar  fields.  Pearson(J 1 4)  extended  their  treatment  by  allowing 
for  anisotropic  velocity  distributions  and  ambipolar  fields.  Baldwin(5) 
considers  wave  propagation  in  inhomogeneous  magnetoplasmas  in  the  low 
temperature  limit,  with  no  resonant  particle  effects.  Analyses  such 
as  these  inherently  lack  generality,  being  based  on  various  approxi¬ 
mations  which  effectively  remove  the  nonlocal  nature  of  the  interactions, 
or  require  low  temperature  or  only  slight  inhomogeneity,  or  rely  on 
quasistatic  conditions.  They  can  be  quite  successful  in  interpreting 
experiments  that  conform  to  the  specialized  assumed  conditions,  although 
the  results  may  often  be  valid  in  only  highly  restricted  ranges  of  the 
parameters,  such  as  near  the  second  cyclotron  harmonic.  Due  recognition 
of  nonlocal  interactions  was  given  by  Buneman(6)in  his  analysis  of  the 
Bennett  pinch.  This  leads  to  an  integral  equation,  rather  than  differ¬ 
ential  equations,  but  one  of  a  rather  intractable  character,  due  to  the 
complexity  of  the  particle  orbits  and  to  the  effects  of  particle  motions 
in  resonance  with  the  disturbances.  The  present  work  deals  with  a  class 
of  systems  describable  by  tractable  integral  equations. 
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II .  OUTLINE  OF  METHOD 


Starting  with  a  specified  particle  distribution  function  , 

for  each  species,  at  time  t  =  0  ,  the  problem  may  be  stated  as  that 

of  determining  its  subsequent  time  development,  f(^c,^,t)  ,  since  this 

function  provides  statistical  information  as  to  the  state  of  the  system. 
The  particles  so  described  are  subjected  to  forces  due  to  the  applied 
magnetic  field,  the  space  charge  forces  of  the  ambipolar  electric  field, 
the  self  magnetic  field  due  to  unbalanced  currents,  and  to  small-signal 
electromagnetic  fields  accompanying  a-weak  disturbance,  either  applied 
externally  or  self-generated.  Collisions  are  excluded  from  consideration 
in  this  paper. 

The  evolution  of  the  initial  distribution  function  is  implicit  in 
the  Boltzmann  equation,  which  states  that  the  density  in  phase  space 
remains  constant  along  the  particle  trajectories,  in  the  absence  of 
collisions.  Regardless  of  how  the  equation  is  formulated,  its  solution 
entails  the  determination  of  the  particle  orbits,  under  the  action  of 
all  the  forces.  The  distribution  at  time  t  is  determined  by  evaluating 
the  initial  one  at  that  point  in  phase  space,  (x  ,v  )  ,  where  a 
particle  arriving  at  x  with  velocity  ^v  at  time  ^  originated  at 
t  =  0  .  Since  the  forces  other  than  the  applied  ones  are  themselves 

determined  from  sources  obtained  by  suitable  averaging  over  the  distri¬ 
bution,  the  overall  problem  is  nonlinear.  This  stage  of  the  analysis 
is  hence  to  be  carried  out  implicitly  to  obtain  expressions  for  the 
particle  orbits,  particularly  the  initial  phase  points,  as  functionals 
of  unknown  force  fields,  to  be  found  subsequently  from  the  functional 
form  of  the  distribution.  Clearly,  it  is  essential  to  arrive  at  a 
fairly  explicit  form  of  the  particle  trajectories  in  order  to  permit 
the  requisite  processing  of  the  distribution  function  to  extract  the 
unknown  fie1ds  consistently. 

Two  intermediate  steps  are  involved  before  the  closing  of  the  set 
of  equations  throrgh  the  extraction  of  self-consistent  fields.  First, 
the  distribution  f(x,v,t)  ,  as  yet  found  only  implicitly  in  terms  of 


unknown  force  fields,  is  to  bo  suitably  weighted  and  averaged  to  obtain 
i  lie  electric  charge  and  current  distributions,  still  In  terms  of  tho 
electromagnetic  fields  of  which  they  are  tho  sourcos.  This  step 
represents  the  determination  of  the  constitutive  parameters  of  the 
medium.  The  next  step  is  the  solution  of  Maxwell’s  equations,  with 
these  charges  and  currents  as  sources,  to  obtain  the  fields  which  exert 
the  originally  assumed  forces.  The  condition  of  self-consistency  then 
expresses  the  dispersion  relation  for  the  medium,  including  the  stability 
of  the  unperturbed  state. 

The  anisotropy  of  the  magnctoplasma  results  in  the  tensor  character 
of  the  constitutive  parameters,  so  that  the  entire  analysis  is  best 
dealt  with  by  matrix  calculus.  This  minor  complication  is  overshadowed 
by  tho  further  one  that,  in  the  inhomogeneous  medium,  the  constitutive 
parameters  are  not  expressible  simply  as  permittivity  or  conductivity 
tensors  but,  more  generally,  as  integral  operators  expressing  the  non¬ 
local  nature  of  the  medium's  sensitivity  to  disturbances.  This  is  a 
direct  result  of  the  transport  of  particles  along  orbits  that  span 
regions  of  varying  properties,  with  different  responses  to  perturbations. 
This  global  sensitivity  ultimately  leads  to  an  integral  equation  for 
Lhe  self-consistent  fields. 

It  is  clear  then  that  for  this  program  not  to  be  frustrated  from 
the  start,  it  is  essential  that  the  particle  transport  be  expressible 
analytically,  and  in  relatively  simple  form,  rather  than  merely 
numerically.  This  is  virtually  impossible  to  achieve  in  a  completely 
self-consistent  manner,  however,  for  the  particle  trajectories  are 
generally  highly  complicated  solutions  to  nonlinear  differential 
equations.  In  particular,  the  orbits  in  the  inhomogeneous  magnetoplasma 
are  not  simply  helical,  because  of  the  action  of  the  ambipolar  electric 
field.  Thus,  the  description  of  even  the  equilibrium  state  is  a 
difficult  task  in  fully  realistic  situations. 

The  key  to  the  analysis  to  follow  is  the  construction  of  a  plasma 
configuration  which  is  characterized  by  simply  expressed  unperturbed 
particle  orbits.  The  class  of  orbits  to  be  admitted  includes  those 
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which  solvo  generalized  harmonic  oscillator  equations.  Although  the 
associated  configurations  aro  found  in  the  general  inhomogonoous ,  two- 
species  enso  not  to  bo  strictly  solf-consistont  and  t into- Invariant ,  it 
will  be  shown  that,  by  proper  choico  of  the  parameters,  the  initial 
state  can  be  made  to  persist  for  relatively  long  pnriods  of  time,  if 
not  indefinitely.  The  effects  which  ultimately  unravol  the  quasi- 
equilibrium  unperturbed  state  can  be  made  second  order  and  act  on  n 
relatively  long  time  scale.  Since  tho  restriction  to  collision-free 
systems  has  already  limited  the  time  duration  of  validity  of  the 
analysis,  there  is  no  inconsistency  in  investigating  the  short-term 
stability  of  a  slowly  varying  unperturbed  state. 

Even  with  simple  analytic  expressions  for  the  particle  orbits,  the 

standard  approach  based  on  the  distribution  function  would  load 

ultimately  to  an  integral  equation  whose  kernel  exhibits  certain 

singularities.  These  aro  associated  with  particles  whose  orbits  are 

in  rcson.incc  with  the  propagating  disturbances  and  arc  essential  to  the 

stability  problem  in  that  they  lead  to  collisionlcss  damping  or  growth 

phenomena.  The  standard  procedure  for  dealing  with  the  singularities, 

(7) 

as  formulated  by  Landau,  involves  complex  contour  integrations  that 
arc  inimical  to  convenient  numerical  analysis. 

To  avoid  the  appearance  of  singular  kernels,  the  method  to  be  used 
describes  the  system  in  terms  of  its  "inverse  phase  space  spectrum," 
i.e.  through  the  Fourier  transform  of  the  particle  distribution  function. 
This  has  a  number  of  additional  advantages,  in  that  it  is  an  algebrai¬ 
cally  simpler  description  which  virtually  eliminates  the  two  intermediate 
steps  of  the  program  outlined  above.  The  extraction  of  the  sources  of 
the  perturbation  fiilds  involves  simply  evaluation  or  differentiation 
of  the  spectrum,  rather  than  integration  of  the  distribution  over 
velocity.  In  addition,  the  transform  formulation  makes  the  solution 
of  Maxwell's  equations  merely  a  matter  of  algebraic  manipulation.  As 
a  result,  the  procedure  can  be  carried  forth  analytically  to  a  final 
integral  equation  with  nonsingular  kernel,  which  may  be  solved 
numerically  in  a  straightforward  manner. 
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As  an  illustration  of  the  method,  tho  integral  oquation  oxpressing 
the  dispersion  relation  of  the  medium  is  obtained  for  a  plasma  with  a 
Maxwellian  distribution  in  velocity  and  a  Gaussian  distribution  in 
space . 
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III.  EVOLUTION  OF  INITIAL  STATE 


For  compactness  of  notation  in  the  initial  manipulations,  phase 
space  will  bo  considered  as  that  spanned  by  the  six-vector  0  , 

comprised  of  the  two  position  vectors  in  configuration  and  velocity 
space,  x , v  .  The  first  task  is  the  determination  of  the  particle 
distribution  function  f(0,t)  ,  given  the  initial  one,  f(0,O)  =  fQ(0)  • 

In  the  absence  of  collisions,  tho  Boltzmann  equation  prescribes 

thnt 


f(0,t)  d60  =  fQ(0o)  &  0O  » 

where  0  is  the  initial  phase  pcint  of  a  particle  whose  phase  point 
is  0  at  time  t  .  This  expresses  particle  conservation,  or  continuity 
in  phase  space  along  tho  particle  tr  ijectorics. 

Consider  any  system  for  which  tho  zero-order  orbit  equation  is 
linear  in  the  phaso;  i.c.  has  *hc  form  of  a  generalized  harmonic  oscil¬ 
lator  equation: 


d0/dt  =  Y  0(t) 


(2) 


It  is  assumed  hero  fiat  the  6  X  6  matrix  Y  is  constant,  determined 

by  both  the  externally  applied  force  field  and  the  internal  space  charge 

field  i  associated  with  the  inhomogeneous  distribution.  It  will  be 

shown  later  that  this  assumption  is  oxset  for  a  homogeneous  magneto- 

plasma  and  very  nearly  exact  for  an  i-*homogeneous  one. 

Under  tho  action  of  the  perturbation,  the  orbit  equation  includes 

an  enharmonic  term,  due  to  the  additional  acceleration  bu.lt  into  the 

six-vector  .  The  trajectory  is  then  specified  by  tho 
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dynamic  equation, 


d0(X)/dX  =  Y  0<X)  +  a  (0(X) , X) ,  0<t)  =  0  ,  (3) 

b 

in  which  the  "initial"  condition  directs  the  orbit  to  the  phase  point 
0  at  time  t  .  What  is  required  in  (1)  is  the  initial  phase  0q  =  0(0). 
An  equivalent  integral  equation  for  the  orbit  is 

0<X)  =  e'Y(t_X)0  -  f  e“Y(T_X)ac(0(T),T)  dT  ,  (4) 

A  6 


and  the  initial  phase  is  given  in  terms  of  its  solution  0(X)  by 


0 


o 


e"YTag(0(T) ,T) 


dT 


(5) 


Specializing  these  exact,  general  equations  to  the  case  of  weak 
perturbations  of  the  harmonic  orbits,  the  initial  phase  0q  =  0Q(0,t) 
is  given  to  first  order  in  the  perturbation  ag(0,t)  by 


0 


O 


YX 

e 


0,t-X)  dX 


(6) 


The  assumption  of  a  harmonic  system  has  led  to  a  simple,  explicit 
expression  for  the  initial  phase  point  of  the  perturbed  orbits  for  any 
weak  disturbance. 

To  complete  the  determination  of  the  time  development  of  the 
initial  distribution,  there  remains  to  relate  the  initial  and  final 
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elements  of  phase  space,  d  0q  and  d  0  .  The  relation  is  given  by 
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mm 


the  Jacobian  of  the  transformation: 


=  |det  90  /S 0 1 


d60 


(7) 


But  the  acceleration  affecting  the  particles  is  quite  generally  solenoidal 
in  velocity,  with  the  consequence  that  the  force  fields  are  such  that 
Tr  Y  =  0  and  (9/90) • ag(0, t)  =0  ,  at  least  in  the  absence  of  colli¬ 

sions.  As  a  result,  (6)  leads,  at  least  to  first  order,  to 

6  6 

det  (90q/90)  =1  ,  so  that  d  0q  =  d  0  .  The  element  of  phase  space 

is  invariant. 

The  time  development  of  the  initial  distribution  is  thus  expressed 
to  first  order  in  the  perturbation  by 

f  (0,  t)  =  fQ^e  Yt  "  J  e^a6^e  •  <8) 

Assuming  further  that  the  initial  distribution  is  an  equilibrium 

one  in  the  unperturbed  system  adds  the  condition  that,  with  ag  =  0  , 

-Yt 

f (0 , t)  =  fQ(e  0)  must  be  independent  of  t  ;  i.e., 


f  (o“Yt0)  =  f  (0) 
o  o 


(9) 


Since  this  is  a  condition  on  the  functional  form  of  the  initial 
distribution  and  holds  for  all  0  ,  (8)  is  thereby  simplified  to 


f (0i t) 


YX  ,  -Y\  .  . 
e  afl(e  0,t-A) 
b 


dX^ 


(10) 
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This  explicit  expression  for  the  time  development  of  an  equilibrium 
distribution  under  a  perturbation  ag(0,t)  states  that  the  distribution 
at  phase  0  at  time  t  is  given  to  first  order  by  evaiuating  the 
equilibrium  distribution  fQ(0)  at  the  slightlv  displaced  phase  0-0  , 

where 


0=0  (0j  t) 


-f 


Y\  ,  ■ 
e  a6(e 


Y\ 


0, t-\)  dX 


(ID 


accumulates  the  perturbations  suffered  along  the  trajectory. 
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IV.  INVERSE  PHASE  SPACE  SPECTRA 


Rather  than  proceeding  to  extract  the  velocity  moments  of  the 
distribution  function  by  integration  of  (10)  to  obtain  the  electro¬ 
magnetic  sources,  the  entire  system  description  will  be  relegated  to 
inverse  phase  space.  This  facilitates  the  extraction  of  velocity 
moments,  and  thence  the  electromagnetic  fields,  and  averts  the  appear¬ 
ance  of  singularities  corresponding  to  resonant  particles. 

Inverse  phase  space  is  spanned  by  the  six-vector  0  ,  composed  of 

the  wave  vector  It  and  the  inverse  velocity  vector  .  The  inverse 
phase  space  spectrum  of  the  distribution  is  just  its  six-dimensional 
Fourier  transform: 


F(0,t) 


■/* 


= 


eie'*d6* 


(12) 


The  spectral  description  is  algebraically  simpler,  replaces  convolutions 

by  products  and  drifts  by  phase  factors,  and  yields  the  velocity  moments 

of  the  distribution  by  merely  evaluating  the  spectrum  and  its  derivatives 

(8) 

at  the  origin  in  inverse  velocity  space.  The  distribution  function 
description  is  best  dispensed  with  entirely;  the  initial  state  of  the 
system  is  to  be  specified  directly  in  inverse  phase  space,  by  F  (0) 

The  evolution  of  the  spectrum  follows  from  that  of  the  distribution 
function.  Since 


f  (0,  t)  =  fQ(0  -  0')  =Jl 


-i0  -(0  -  0  ) 

FO(0O)  e  °  d60o/(2*)6  ,  (13) 
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expansion  to  first  order  In  0  and  Fourier  transformation  yields  foi 
the  perturbed  spectrum 


F<e 


,t)  =  Fo<6)  +  i/Fo(eo)0o-e'(0-eo,t)  d6©^*)6  ,  U4) 


where  0  (0,t)  is  the  Fourier  transform  of  0' (0,t)  .  From  (11), 

this  is  given  by 


0*  (0.  t)  =  f  eYXfai6(e~YX0,  t-X)  e10'0  d60  dX  , 


(15) 


which  indicates  that  the  relevant  perturbation  quantity  is  the  acceleration 

in  inverse  phase  space,  A  (0, t)  ,  the  Fourier  transform  of  a_(0,t) 

6  6  ® 

Since  d  0  =  d  (e  0)  ,  (IS)  is  expressible  as 


0*  (0,t)  =  f  eYX  Afi(0 

»  n 


eYX,  t  -  X)  dX 


(16) 


6  6  YX 

and  a  similar  change  of  variable,  using  d  0Q  =  d  (G„e  )  and  the 


Yt. 


fact  that  the  equilibrium  condition,  (9),  translates  into  F  (0)  =  F  (0e  ), 

o  o 

converts  (14)  into 


F  (0, t)  =  i  f  f F  (0.-0  )(0.-0  ) ’ A  (0  ,t-X)d60  dX/(2«)6  .  (17) 

1  J  Jo\oko6o  o' 
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Here,  F^(0,t)  is  the  perturbation  of  the  spectrum  and 


et  =  0  e 


Yt 


(18) 


may  be  viewed  as  tht  orbit  in  inverse  phase  space. 

With  the  perturbation  of  the  spectrum  now  expressed  as  a  convolution 
in  time,  a  final  Laplace  transformation  is  indicated,  which  will  also 
reduce  Maxwell's  equations  to  algebraic  ones.  Denoting  the  Laplace 
transform  operation  by  G(s)  =  L  g(t)  ,  the  convolution  theorem  yields 
for  the  transformed  spectrum  perturbation, 

F  (0 , s)  =  iL  ft  <0-0  )(0  -0  )-A.(9  .s)d60  / (2n) 6  .  (19) 

1  ■~ej  ototoeo  o 

This  relation  specifies  the  operations  to  be  carried  out  on  the  equilib¬ 
rium  spectrum  F  (0)  to  obtain  the  system  response  to  the  perturbation, 
o 

Fig.  1  traces  the  methods  of  calculating  the  response  of  the  system 
to  the  perturbing  fields,  in  the  phase-time,  inverse  phase-time,  and 
inverse  phase-frequency  domains. 


9 
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V.  SELF-CONSISTENT  FIELDS 


The  acceleration  associated  with  a  perturbation  consists  of  any 
weak  externally  applied  field,  such  as  an  impulse  or  a  signal  wave, 
together  with  the  reaction  provided  by  the  internal  electromagnetic 
fields:  Ag(9,s)  =  Ag  Xt(9,s)  +  A6em(9,s)  .  In  particular,  the  plasma 

response  to  a  weak  impulse  which  imparts  velocity  v  uniformly  to  the 
particles  is  obtainable  by  setting 


ag(0,t)  = 

0  ' 

V 

6(t) 

Ag ( 0 , s )  = 

'0  ‘ 

V 

.  P. 

.  P. 

(2Tt)66(0) 


(20) 


so  that,  from  (19)  and  the  equilibrium  condition, 


F  (0,s)  =  iL  F  (0  e  r)  0 


~s  o 


0 

VPJ 


+  Fiem(0,s) 


iFQ(0)  0 ■ (s  -  Y) 


-1 


L  P. 


_  em.. 

+  (0,s) 


(21) 


There  remains  to  obtain  the  electromagnetic  interaction. 

Three  steps  are  involved  in  expressing  the  internally  generated 


em. 


perturbation  F^'^O.s) 


in  terms  of  the  total  perturbation  F1(0,s) 


By  virtue  of  the  present  transform  formulation,  each  of  these  has  been 
reduced  to  algebraic  manipulation.  The  first  step  is  the  extraction 
of  the  field  sources  from  the  perturbed  spectrum;  Maxwell's  equations 
then  yield  the  r.f.  electromagnetic  fields  generated  by  these  sources. 
Finally,  the  internally  generated  acceleration  is  expressed  in  terms 
of  the  perturbation  through  the  Lorentz  force.  This  closes  the  equation 
self-consistently . 
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If  the  system  Includes  N  particles  of  mass  m  and  charge  (-e) 
distributed  in  phase  space  as  f(#,t)  ,  then  with  the  normalization 

P  (0)  =  1  ,  the  perturbation  charge  and  current  densities  contributed 

by  that  species  are  given  by  their  Fourier-Laplace  transforms  as 


p^k.s)  =  "e  N  ViS>°'s)  ’  £i(iS's)  =  ie  N  9  F1<^»°»s>/SA  •  (22) 


Note  that  particle  conservation  is  expressed  by 


k 

s 


SF^k.O.s) 
1  ~ 

9A 


=  F  (k,0,s) 

1  ~ 


(23) 


Maxwell's  equations  combine  into  a  wave  equation,  which  transforms 


to 


k  X  (k  X  a)  =  S2  a  -  (e/meQ)  s  Jj/k.s) 


(24) 


where  a(k,s)  =  (-e/m)  E(k,s)  is  the  transform  of  the  acceleration 
suffered  by  a  particle  subjected  to  only  the  r.f.  electric  field,  and 
J1(k,s)  is  the  total  current  density  due  to  all  the  species.  Continuing 
to  consider  only  the  contribution  of  the  one  species,  (24)  may  be 
expressed  in  terms  of  the  perturbed  spectrum  by 


c2(k  k  -  k2)-a  =  s2a  -  iof  s  9F  (k,0,s)/9A 

Vw  ^  -V  O  A 


(25) 


where  a>2  =  Ne2/m€  is  the  product  of  the  average  squared  plasma 
frequency  and  the  volume  containing  the  N  particles.  Together  with 
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(23),  this  equation  incorporates  Gauss'  law, 


k-a  = 


ko  F  (k, 0, s) 

O  1  ~ 


(26) 


so  that  the  solution  of  the  wave  equation  is 


Q(k, s) 


io) 


2,2  2 

s  +  k  c 


(°2~  +  8 


s) 


(27) 


By  Faraday's  law,  the  transform  of  the  acceleration  due  to  the  r.f. 
magnetic  field  interaction  is  v  X  (ik  X  a)  ,  so  that  the  total  Lorentz 
acceleration  is  given  in  inverse  phase  space  by 


'o' 

Q 

/k  k\ 

A  (k,A,s)  = 
o  ~ 

1 

(2k) 

a  6(A)  +  (  -  a  -  a  -  J-6 '  (A) 

(28) 


It  should  be  remarked  that  this  formulation  neglects  relativistic 

2 

corrections  of  order  (v/c)  but  retains  all  terms  of  order  (v/c) 

Introducing  A_(k  ,A  ,s)  into  (19)  and  integrating  over  A 
b  ~o  ~o  ~o 

leaves,  after  some  manipulations, 


em 


(0,s)  = 


iL 

~s 


/['  * 


A^  + 


(k 

~o 


/s>] 


a(k  ,s)d3k  /(2tt)3  ,  (29) 


where 


9t  - 


(It,, A,) 


=  0eYt  =  (k,A)eYt 


Ft  = 


F  (k.-k  , A  ) 
o  ~t  ~o  ~t 


(30) 
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and 


G  (k,A)  =  A  x 

/N^O  ^  ~ 


y.'  (k,A) 

O  ~ 


G  =  G  (k-k, A, ) 
~o  ~t  ~o  ~t 


Th.  first  and  second  terns  In  the  Integral  In  (29)  represent  the  r.f. 
electric  and  magnetic  Interactions,  respectively,  In  the  plasma.  The 
equation  gives  the  spectrum  perturbation  accompanying  the  electric 
field  represented  by  2(k.s>  ,  while  (27)  yields  the  field  associated 

with  the  perturbation  In  the  spectrum.  In  the  absence  of  external 
excitation,  the  self-const-tent  fields  are  the  solution  of  the  homo¬ 
geneous  integral  equation  formed  by  combining  (27)  and  (29),  with 
_  em  _  p  impulse  response  Is  given  by  the  solution  of  the 

Inhomogeneous  Integral  equation  obtained  by  combining  (27),  (29)  and 

(21) 

The  primary  objective  of  reducing  the  problem  to  th.  solution  of  a 
nonsingular  Integral  equation  has  thus  been  attained.  The  kernel  o 
the  final  integral  equation  Is  essentially  the  coefficient  of  2«o-> 
in  the  integral  In  (29):  It  Is  formed  from  the  function  F0(k,~ 
representing  the  equilibrium  spectrum  by  explicit  operations  specl  led 
in  (30) ,  (31) ,  and  (29) .  It  Is  generally  a  well-behaved,  nonslngu  ar 
kernel.  The  Integral  equation  ostensibly  determines  the  plasm,  frequency, 
„  ,  as  Its  eigenvalue.  Since  this  Is  actually  a  specified  parameter 

of  the  system,  however,  the  determination  of  the  eigenvalue 
expresses  the  self-consistency  condition,  or  the  dispersion  re  at Ion 
restricting  s  and  k  .  The  fact  that  %  Is  a  real,  pcs  ve 
quantity  Imposes  restrictions  on  the  complex  frequency  s  ,  which 
stipulates  the  stability  of  the  equilibrium  state. 

Summarizing  the  procedure  for  determining  the  stability,  dispersion 
relation,  or  impulse  response  of  a  collection  of  particles,  the  sys  em 
is  first  described  by  some  Inverse  phase  space  spectrum  Fq(9)  , 

each  species.  The  force  fields,  Including  those  externally  app 
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and  the  internal  ones  which  maintain  the  distribution,  are  described  in 
terms  of  the  accelerations  produced  and  are  assumed  to  compose  a  gener¬ 
alized  harmonic  oscillator,  specified  by  the  constant  6X6  matrix  Y 
which  validates  (2) .  The  restrictions  on  these  initial  specifications 

are  the  solenoidal  field  condition  implied  by  Tr  Y  =  0  and  the  equi- 

Yt 

librium  condition  expressed  by  F  (0e  )  =  F  (0) 

o  o 

The  vector  G  (0)  =  A  X  SF  /3A  determines  the  magnetic  interaction 
~o  ~  o  ~ 

and  essentially  measures  the  anisotropy  of  the  equilibrium  velocity  dis- 

Yt 

tribution.  The  orbit  in  inverse  phase  space,  (lc .  , A .  )  =  (k,A)e  ,  may 

Yt~  ~ 

be  evaluated  by  standard  matrix  methods;  the  relation  e  =  L  (s-Y) 
probably  yields  this  most  easily.  Then  the  perturbation  of  the  spectrum 
is  given  in  terms  of  the  first-order  electric  field  by  (29),  while  the 
latter  is  given  in  terms  of  the  former  by  (27).  The  combination  forms 
an  integral  equation  for  a(k,s)  ,  whose  eigenvalues  determine  the 
dispersion  relation  among  s,  £  ,  and  the  plasma  frequency,  o)Q  .  The 

response  to  applied  signals  is  obtained  from  an  inhomogeneous  integral 
equation,  formed  as  in  (21)  for  the  impulse  response. 

Fig.  2  depicts  the  processing  of  the  given  equilibrium  spectrum 
and  acceleration  matrix  required  to  arrive  at  the  equations  relating 
the  spectrum  and  field  perturbations. 

To  account  for  perturbations  of  several  species,  the  individual 
fields  of  each  may  be  superposed  to  form  a  .  Any  small  deviations 
from  the  harmonic  oscillator  equation,  (2),  may  be  treated  as  a  further 
perturbation,  as  in  (3) .  „ 
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VI .  GAUSSI AN-MAXWELLIAN  COLUMN 


Before  llustrating  the  method  of  sotting  up  the  integral  equation 
with  an  application  to  a  particular  inhomogeneous  magnetoplnsma,  it  may 
be  noted  that  for  the  homogeneous  case,  the  integral  equation  is  elim¬ 
inated,  reducing  to  an  algebraic  relation.  This  reduction  is  made 
possible  physically  by  the  fact  that  the  particles  then  traverse  only 
regions  which  respond  to  perturbations  just  as  does  the  one  in  the 
immediate  vicinity  of  the  initial  position,  so  that  the  global  sensi¬ 
tivity  is  reduced  to  a  local  interaction.  Formally,  a  homogeneous 
system  is  described  by  a  spectrum  of  the  form 


F  (6)  =  F  Ck,A)  =  F(A)  (2jr)3  6(k) 

O  O  '■v  '■v 


(32) 


whereupon  (29)  reduces  to 


em,  a 

F.  (0,sj 


iL 

~s 


[F<i, 


A  -a(k  ,s) 


+  G(A  ) 


‘r  x  “<£,>■>] 


(33) 


where  0(A)  =  A  X  dF/dA  .  The  self-consistent  field  equation  is  then 
obtained  by  combining  this  with  (27),  algebraically. 

A  fairly  realistic  example  of  a  warm,  inhomogeneous  magnetoplasma 
is  provided  by  one  described  by  a  Gaussian-Maxwellian  spectrum: 

Fo(6)  =  exp  (-  j  8  R  0)  .  (34) 


Here,  R  is  a  constant,  symmetric  6X6  matrix,  which  partitions  into 
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four  3X3  submatrices  as 


R  = 


(35) 


incorporating  the  parameters  that  specify  the  mean  spatial  extent  of 
the  plasma,  its  thermal  velocity,  and  its  drift  motion.  As  a  matter  of 
interest,  for  ease  of  visualization  but  superfluous  to  the  calculation, 
this  spectrum  corresponds  to  the  particle  distribution  function 


fo<*> 


exp  (-  ^  0R  1  0)/( 2jt)3 


y2 

(det  R) 


(36) 


which  represents  a  collection  of  particles  with  a  Gaussian  spatial 
distribution  and  a  Maxwellian  velocity  distribution.  The  particle 
drift  is  linear  in  position,  u(x)  =  R„  R~*  x  ,  and  the  temperature 

~  ~  VX  XX  ^  2 

is  uniform,  specified  by  the  mean  squared  thermal  velocity,  T/m  =  v^  : 


2  1  „  /n 
v.  =  -  Tr(R 
0  3  vv 


R  R"1  R  ) 

vx  XX  XV 


(37) 


These,  and  other  velocity  and  spatial  moments  of  the  distribution,  are 
z  vidily  obtainable  by  differentiation  of  the  spectrum  and  evaluation  at 
the  origin  of  A-  or  k-space. 

The  equilibrium  force  field  acting  on  the  plasma  is  assumed  to  be 
expressible  by  the  constant  6X6  matrix,  Y  ,  relating  d0/dt  to  0  , 
forming  a  generalized  harmonic  oscillator.  This  partitions  into  four 
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3X3  submatrices  as 


Y  = 


0 

Y 

vx 


(38) 


The  submatrix  Yyy  is  to  have  zero  trace  to  conserve  the  element  of 
phase  space.  It  represents  an  acceleration  linear  in  velocity,  as  is 
the  case  for  a  magnetoplasma.  Space  charge  fields  are  to  be  accounted 
for  through  Yyx  ,  which  admits  an  acceleration  linear  in  position, 
representing  a  particular  class  of  ambipolar  fields. 

The  condition  for  the  unperturbed  distribution  to  be  an  equilibrium 
state  is  here  expressed  by 


Fo<eeYt>  =  exp  (-  |  0eYt  ReY  t  0)  =  Fq(0)  =  exp  (-  |  0  R  0)  ,(39) 


where  the  prime  denotes  the  transpose  of  the  matrix.  By  inspection  of 
the  quadratic  form,  it  is  readily  seen  that  the  equilibrium  condition 
reduces  to  the  requirement  that  YR  be  antisymmetric. 

Specializing  to  a  uniaxially  symmetric  magnetoplasma  column,  uniform 
in  the  axial  direction  but  transversely  Gaussian,  reduces  the  6X6 
matrices  to  4  X  4,  eliminates  the  integration  over  the  axial  wave 
number,  and  requires  N  to  be  redefined  as  the  axial  linear  particle 
density.  A  suitable  combination  of  R  and  Y  matrices  that  satisfy 
all  the  conditions  is  specified  for  electrons  by  the  2X2  submatrices 


2  „ 

R  =  a  J 
xx  e 


vv 


=  V2  I 


XV 


R  =  0 
vx 


t 


Y 

vx 


2 


“CD 


e 


I 


Y 

vv 


CO  x 

ce 


X 


0  -1 

1  0 


(40) 
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The  last  of  these  effects  the  cross  product  operation  with  the  axial 
magnetic  field;  coce  is  the  electron  cyclotron  frequency.  There  is 
thereby  described  a  rragnetoplasma  column  with  Gaussian  radial  density 
profile  of  effective  radius  ag  and  Maxwellian  velocity  distribution 
with  thermal  velocity  Vg  .  Equilibrium  obtains  when  o>  =  v  /a  • 
this  transit  frequency  corresponds  to  a  time  scale  of  the^rder  of  the 

time  it  would  take  a  thermal  electron  to  traverse  the  effective  radius 
of  the  column. 

Fig.  3  presents  the  essential  characteristics  of  the  Gauss ian- 
Maxwellian  magnetoplasma  column  here  considered. 

The  configuration  defined  by  (40)  Includes  no  particle  drift  but 

implies  an  amblpolar  field  providing  a  radially  Inward  acceleration 

(  toe  r)  which  maintains  the  Gaussian  profile.  This  field  Is  to  he 

provided  by  the  Interaction  with  the  Ion  distribution.  But  this  radially 

increasing  amblpolar  field  which  is  to  make  the  equilibrium  configuration 

.  harmonic  oscillator  strictly  requires  charges  at  Infinity,  or  violates 

overall  neutrality.  This  system  Is  hence  not  strictly  resile. ble  on  a 

permanent  basis  of  macroscopic  self-consistency.  It  Is  however  possible 

to  maintain  the  desired  conditions,  Including  chsrge  neutrality,  at 

least  temporarily,  by  an  appropriate  selection  of  the  Ion  distribution 
as  follows. 

If  the  same  number,  N  ,  of  ions  per  unit  length  be  distributed 
with  a  Gaussian  profile  of  effective  radius  ^  slightly  smaller  than 

3e  ’  then  the  res>ll‘lng  amblpolar  electric  field  would  provide  an 
inward  acceleration  to  the  electrons,  of  magnitude 


2jrr  |exp  <-  2  ~2>  "exP  (-  o  -7?) 
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nocif iciTiocs.  The  correct ions  arc  the  s*all*uv  the  WsMv  wwUvM't  the 
radial  distribution  and  the  lower  The  tewjwatuve* 

The  asibipolar  field  which  he  re  restrains  t  he  elec  t  VOW*  tends  1  e 
disrupt  the  ion  distribution,  which  hence  fails  to  center**  te  the 
postulated  harmonic  oscillator  dynamics.  The  initial  slat*'  here 
described  is  hence  only  a  quasi-equilibrium  state  tut  the  lurin'  ton 
to-electron  mass  ratio  ensures  that  the  time  dove l opmeu t  of  the  *tts 
tribution  can  be  quite  slov.1  compared  to  the  time  seal*'  of  t tie  i\f> 
perturbations  here  considered.  Since  collisions  were  excluded  nt  III*' 
outset,  the  validity  of  the  formulation  extends  In  time  for  only  a 
mean  collision  time.  It  follows  that  the  Inch  or  permanent  equilibrium 
is  no  severe  restriction  here  since  all  that  need  be  Imposed  Is  Hull 
the  parameters  allow  the  Gaussian  distribution  In  persist  ensoul  I tt 1 1 V 
undistorted,  in  the  absence  of  the  port  itrlm!  tons ,  for  one  noil  la  Ion 
mean  free  time.  The  description  Is  exact,  ioc  \\  pi  mm  him  wltlelt  la  ell  her 
cold  or  homogeneous;  it  is  a  close  approx i mat lou  otherwise)  provided 
the  inhomogeneity  and  temperature  uro  not  extreme  and  llu*  muss  ratio  Is 
large. 


With  the  understanding  that  the  system  is  being  observed  only 
during  that  time  interval  in  which  it  behaves  as  a  harmonic  oscillator 

in  equilibrium,  the  analysis  proceeds  by  specializing  the  general  result 
to  the  spectrum 


F  (6)  =  exp  -  -  (k2af  +  A2v2)2n  6(k  ) 

°  £  l  e  e  II 


(43) 


and  to  the  orbit  in  inverse  phase  space  given  axially  by 


(k  ,A  )  =  (k  ,A  +  k  t) 

t  t  II  II  II  II 


(44) 


corresponding  to  free  particle  motion  along  the  magnetic  field,  and 
transversely  by 


<kf\\  =41<k'A) 


-  -1 


-I 


03_ 


S-fflceX 


(45) 


which  involves  the  magnetic  and  ambipolar  fields.  Since  the  velocity 
distribution  is  here  isotropic,  G  (9)  =  0  and  there  is  no  first- 
order  r  f.  magnetic  interaction.  The  spectrum  perturbation  is  hence 
given  in  terms  of  the  perturbing  acceleration  a(k,s)  by 


F^O.s) 


2  AtVe 


a 


2  is  6 


JL-y(K 

~t  ~  ~t 


(46) 
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where  ^(k)  is  a  Gaussian  transform  of  Q(k)  : 


Y(k,s) 


exp  -  i  (k-k  )2a2 
*  «s.>  *jo  e 


a(k  , 1 

J  ~  ~o 


,s)d2k 


(47) 


The  axial  homogeneity  gives  the  axial  component  k  the  role  of  a 
parameter.  Note  that  in  the  fully  homogeneous  case,  a  ”*  <»  and 
Y(k)  =  a(k)  .  6 

~  ~  «v  <v 

The  integral  equation  for  the  self-consistent  field  a(k,s)  is 
now  obtained  by  combining  (46)  with  (27).  Defining  the  plasma  frequency 
cup  in  terms  of  the  axi£l  density  by 


2  .  2 
=  Ne  /2jta  me 
e  o 


(48) 


yields  the  result 


a(U,s) 


0) 


2 

s  + 


.2  2 
k  c 


(°  ~ + s  4) 


L 

~s 


1  .2  2 

■  2  V. 


=  0.(49) 


A=0 


A  more  explicit  form  of  this  integral  equation  is  derivable  by  use  of 
the  equilibrium  condition,  which  here  imposes 


Substitution  in  (49)  reduces  the  integral  equation  to 
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An  alternate  form  involving  an  exponential  rather  than  a  Gaussian  kernel 
follows  from  this  by  defining 


£05’ s> 
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(52) 


The  unknown  is  now  the  convolution  of  the  electric  field  with  the  density 
distribution  and  the  integral  equation  becomes 
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Further  insight  into  the  properties  of  these  integral  equations 
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may  be  gained  by  examining  the  explicit  orbit  in  inverse  phase  space 
Developing  (45)  yields  for  the  transverse  motion  complementing  (44), 

ktx  =  kx  Z(t)  +  \  dZ/dt  •  Atx  =  <-l/o£)dkt  /dt  ,  i 

where  Z(t)  is  determined  by  a  matrix  partial  fraction  expansion  of 
(45)  as 


“2 


L  Z(t)  = 

~S  CD,  +  CD- 


-1 

(s  +  CD^X)  + 


“l  -1 

- ; -  (s  -  cd„X) 

03,  +  CD„  2 


(55) 


or 


Z(t)  = 


“2 

<-“>!*»  -  5^—  «p  stX)  ■ 


(56) 


This  combines  two  rotations,  in  opposite  senses,  at  rates  whose  difference 
is  the  cyclotron  frequency  and  whose  geometric  mean  is  the  transit 
frequency : 
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1  ce 
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(57) 


The  trigonometric  form  of  (5G)  is 
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which  indicates  that  the  kernel  in  (53)  actually  involves  Bessel 

(9) 

functions.  At  A  =  0  ,  for  k*k  =  kk  cos  A  , 
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Introducing  the  orbit  into  (53)  reduces  the  integral  equation  to 
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is  merely  linear  in  k 

/vO 
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The  above  analysis  illustrates  the  method  of  arriving  at  the 
integral  equation  governing  an  inhomogeneous  magnetoplasma  and  the  forms 
which  that  equation  can  take.  The  numerical  analysis  which  remains  to 
be  carried  out  on  the  equation  to  extract  the  dispersion  relation  from 
the  eigenvalues  involves  computations  on  only  well-behaved  kernels. 

While  the  final  numerical  processing  of  the  equations  to  reveal  the 
wave  supporting  and  stability  properties  of  the  plasma  here  considered 
is  beyond  the  scope  of  this  expository  work,  a  few  general  remarks  on 
the  character  of  the  integral  equations  may  be  in  order. 

Several  distinct  but  equivalent  forms  of  the  final  integral 
equation  have  been  presented.  The  kernel  may  be  Gaussian,  as  in  (51), 
or  essentially  exponential,  as  in  (60),  or  a  superposition  of  Bessel  ’ 
functions,  as  when  the  expansion  (59)  is  substituted  into  (60).  The 
choice  permitted  by  the  commutativity  of  the  integral  operator  and 
the  Laplace  transformation  in  any  of  these  forms  provides  further 
latitude  in  the  selection  of  the  equation  to  be  analyzed  numerically. 

For  example,  in  (60),  the  Laplace  transformation  operates  on  the  time 
function  indicated,  which  involves  H^.s.t),  whose  Laplace  transform 
is  trivially  deduced  from  (55) ,  modified  by  what  is  purely  an  exponential 
function  of  t  if  (59)  is  introduced  into  (60).  ;[n  the  absence  of 
the  factor  exp  (-  |  k^t2)  >  th®se  exponentials  would  result  in 
merely  a  shift  in  the  frequency  domain,  with  the  Laplace  transform  of 
H  evaluated  at  [s  -  i(n  ^  +  m  o>2)j  .  The  effect  of  the  exponential 

involving  the  axial  wave  vector  k |  is  to  introduce  a  Gaussian  factor 
in  the  transformation,  leading  to  the  familiar  error  function,  or  plasma 
dispersion  function,^  ^  and  to  Landau  damping. 

2  2he  connection  to  the  physical  system  is  T>ade  through  the  factor 
(“p/^e*  in  (6°)*  which  plays  the  role  of  the  eigenvalue  of  the  integral 
equation  but  effectively  relates  the  parameters  s  and  k  ,  thereby 
expressing  the  dispersion  relation.  The  stability  of  the  Equilibrium 
state  is  determined  by  those  locations  of  s  in  the  complex  plane 
which  permit  the  eigenvalue  to  be  real  and  positive. 
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The  two-dimensional  character  of  the  Integral  equations  can  be 
overcome  in  this  case  by  taking  advantage  of  the  symmetry  of  the  column 
to  eliminate  the  azimuthal  integration.  By  use  of  (59) ,  this  reduces 
the  double  summation  to  a  single  one,  as  well  as  leaving  only  a  one¬ 
dimensional  integration  over  the  radial  wave  number  k  ,  suitable  for 

o 

digital  computation. 

Various  approximations  of  both  a  mathematical  and  physical  nature 

may  be  imposed  on  the  equations.  The  common  quasistatic  approximation 

may  be  introduced  by  formally  letting  c-*  a>  In  (60),  this  reduces 
2  2  2  2 

H0c,kQ, s , t)/(s  +  k  c  )  to  (k  Jt/k  )dZ/dt  and  considerably  simplifies 

the  equation.  Other  simplif ications  are  effected  by  selecting  certain 
directions  of  wave  propagation  and  polarization  with  respect  to  the 

magnetic  field,  particularly  k  =  0  .  The  combined  case  is  treated 

li 

more  thoroughly  in  the  Appendix. 

For  a  mildly  inhomogeneous  and  warm  plasma,  the  transit  frequency 

to  is  low  and  the  effective  radius  a  is  large.  Suitable  Taylor 
©  © 

expansions  of  the  kernels  then  permit  the  determination  of  at  least  the 

first-order  deviations  from  the  homogeneous  case.  It  is  already  clear 

from  (55-59)  that  to  the  very  lowest  order,  the  cyclotron  frequency 

and  its  harmonics  are  replaced  by  the  two  frequencies  to.  ,cn_  and  all 

their  combination  frequencies;  cc  differs  from  u)  by  cc.  ,  which 

^  C©  X 

is  small  for  low  transit  frequencies.  Beyond  that,  approximate  results 

are  obtainable  by  noting  that  the  Gaussian  factor  in  (47)  becomes 

sharply  peaked  as  a  °°  .  A  Taylor  expansion  of  a(k  )  about 

e  ~  ~o 

k  =  k  yields 


y(k)  =  a(k)  + 

2a 

e 


(62) 


When  this  is  introduced  into  (51),  there  results  a  differential  equation 
as  an  approximation  to  the  integral  equation. 

Finally,  it  may  be  surmised  from  a  qualitative  examination  of  the 
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integral  equations  with  Gaussian  kernels  that  whereas  plane  wave  solutions 
of  the  form  6(k  -  kg)  are  appropriate  in  the  limiting  infinite, 
homogeneous  case,  the  solutions  for  a(k)  will  in  the  general  case 
involve  a  spectrum  of  k  vectors,  sharply  peaked  about  some  signal 
wave  vector  kg  The  dispersion  relation  then  relates  the  signal 
frequency  to  this  central  vector  of  the  wave  packet,  which  is  however 
comprised  of  a  broadened  spectrum  of  k  vectors. 
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VII.  CONCLUSIONS 


The  self-consistent  field  problem  for  a  cl&ss  of  anisotropic, 
inhomogeneous  particle  distributions  has  been  reduced  to  that  of  solving 
an  integral  equation  with  a  well-behaved  kernel.  Tractability  has  been 
achieved  mainly  by  considering  a  class  of  systems  for\which  the  equili¬ 
brium  particle  orbits  can  be  simply  described  analytically.  In  partic¬ 
ular,  systems  whose  dynamics  correspond  to  small  deviations  from 
generalized  harmonic  oscillation  have  been  analyzed.  The  particle 
propagator  then  has  the  form  exp  Yt  ,  making  it  relatively  simple  to 
evaluate  and  accumulate  the  perturbations  along  the  orbit. 

The  approach  can  evidently  be  generalized  to  other  systems  with 
known  and  simply  expressed  unperturbed  particle  orbits.  A  rather 
direct  generalization  might  be  to  certain  time-varying  systems,  in  which 
the  propagator  has  the  form  exp  £  Y(t)  dr  ,  provided  the  acceleration 
matrix  Y(t)  commutes  with  its  integral.  A  more  sweeping  generalization 
might  be  made  to  orbit  equations  quadratic  in  the  phase,  rather  than 
merely  linear.  The  orbit  then  satisfies  a  Riccati  equation  and  the 
systems  so  described  admit  ambipolar  fields  more  general  than  those 
considered  herein. 

Simple  particle  orbits  permit,  simple  recipes  for  the  evolution  of 
an  arbitrary  initial  particle  distribution  function  in  response  to 
arbitrary  perturbations.  The  relevant  physical  quantities  are,  however, 
various  moments  of  the  distribution,  or  other  expectation  values.  Their 
extraction  by  integration  over  velocity  space  is  not  only  cumbersome 
but  complicated  by  classes  of  particles  whose  orbits  are  in  resonance 
with  the  propagation  of  disturbances.  These  difficulties  are  avoided 
by  describing  the  system  at  the  outset  by  its  spectrum  in  inverse  phase 
space.  This  spectrum  is  the  characteristic  function  of  the  distribution, 
or  the  expectation  value  of  exp  i 0-0  ,  and  exhibits  all  the  moments 

by  its  behavior  in  the  immediate  vicinity  of  the  origin  of  inverse  phase 
space. 
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The  perturbation  of  the  spectrum  is  expressed  in  terms  of  the 
Four ier-Lapl ace  transform  of  the  perturbation,  which  is  algebraically 
related  to  the  sources  of  the  electromagnetic  fields,  which  are  in 
turn  given  directly  by  the  perturbed  spectrum.  The  system  of  equations 
is  thus  closed  self-consistently  through  an  integral  equation  whose 
kernel  is  essentially  the  unperturbed  spectrum  evaluated  along  the 
orbit  in  inverse  phase  space.  The  equation  admits  arbitrary  combinations 
of  equilibrium  spectra  and  the  associated  harmonic  oscillator  force 
fields  and  involves  well-behaved  kernels , well-suited  for  final  numerical 
analysis 

The  case  of  a  Gauss ian-Maxwe Ilian  plasma  has  been  treated  in  greater 
detail,  by  way  of  illustration.  The  plasma  configurations  conforming 
to  the  exp  (-  |  0  R  8)  spectrum  have  not  been  exhausted  herein.  One 
which  requires  no  ambipolar  electric  field  to  maintain  its  density 
profile  is  a  column  described  transversely  by 


R 


4 


coco 

c 


(63) 


This  has  a  Gaussian  radial  profile  with  effective  radius  a  ,  a 

Maxwellian  velocity  distribution  of  temperature  T  =  m  a2co(co  -  co)  . 

c 

and  an  azimuthal  drift  corresponding  to  rotation  at  the  rate  co  about 

the  axis  of  the  column.  Except  for  diamagnetic  effects  of  the 

azimuthal  drift  current,  this  column  is  in  equilibrium  in  an  axial 

magnetic  field  specified  by  the  cyclotron  frequency  co  .  The  foregoing 

c 

analysis  applies  to  this  configuration,  with  R  replaced  by  (63), 

Yvx  =  0  ,  Yyv  =  cocX  .  It  may  also  be  noted  here  that  displacements 
through  x  and  drifts  at  velocity  v  can  be  simply  incorporated  into 
the  inverse  phase  space  spectra,  through  a  phase  factor  exp  i0-0 
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The  illustrative  example  presented  here  affords  an  opportunity  to 
formulate  integral  equations  for  a  plasma  in  which  temperature,  anisotropy, 
inhomogeneity,  particle  transport,  and  ambipolar  fields  are  all  to  be 
accounted  for.  There  appears  to  be  no  insurmountable  obstacle,  other 
than  increased  complexity,  to  generalizing  the  formalism  to  include 
collisions  and  relativistic  effects. 
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APPENDIX 


TRANSVERSE  OSCILLATIONS  OP  MAGNETOPLASMA  COLUMN 


Under  the  quasistatic  approximation,  the  special  case  of  a  self- 

consistent,  azimuthally  symmetric  transverse  field  pattern,  which  may 

represent  the  cutoff  condition  for  axial  propagation  along  the  Gaussian- 

Maxwellian  magnetoplasma  column,  is  sufficiently  simple  to  analyze 

quite  thoroughly.  With  reference  to  (52,60,61),  the  conditions  are 

expressed  by  c  -•  «  ,  k  =  0  ,  a(k,s)  =  nh(u,s)  ,  where  n  is  a  unit 
r  "  ||  rSm>  ^  *** 
vector  along  k  and  u  =  ka  .  It  then  suffices  to  set  A  =  0  ,  so 

~  e 

that  the  orbit  is  given  by  k  =  kZ(t)  and  the  integral  equation  reduces 
to 
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Choosing  polar  coordinates  with  polar  axis  along  ri  and  azimuth  defined 

by  n  =  exp  (rtX)n  ,  the  azimuthal  integration  becomes,  using  (56), 

M>  F  ~ 


uu  nZn 
o~  ~o 


e 


d0 


exp (uu 

o 


cos(0  -  03^t) 


+ 


“2 

—  cos(0  +  03_t) 

0)3  2 


d^ 


=  i  (un(t)u  ) 
o  o 


(A-2) 
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Consequently,  the  Integral  equation  collapses  to 
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2  2 

where  X  =  0)^/0^  is  the  eigenvalue. 

By  use  of  the  multiplication  theorem  for  Bessel  functions/12*  it 
is  possible  to  reduce  the  kernel  of  this  equation  to  degenerate  form. 
Since 
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the  equation  is  equivalent  to 
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The  degeneracy  now  permits  the  reduction  of  the  integral  equation  to  an 
equivalent  infinite  set  of  simultaneous  equations  of  the  form 


H  <s)  =  \  S  G  (s)  H  (s)  ,  (A-8) 

n  nm  m 
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with 


G  (s)  = 
nm 


L  e  V"*"  ks  Tt  - 
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.  (A- 9) 


These  matrix  elements  may  be  evaluated  explicitly,  by  virtue  of  the 
relations 
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and,  in  terms  of  Laguerre  polynomials, 
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and  the  resonance  condition  is 


det  (I  -  \G)  =  0  , 


(A-13) 
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which  relates  the  complex  frequency  s  to  the  plasma  frequency  co 

22  22  l/2  ^ 

in  \  =  05/05  .  Clearly,  the  harmonics  of  05o  =  (to  +  4o>  )  ' 

p'  e  3  ce  e 

replace  the  cyclotron  harmonics  as  the  significant  frequencies. 
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In  (p, t)  ,  (0,t)  ,  and  (0,s)  spaces. 
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EVOLUTION  OF  INTEGRAL  EQUATION 


Fig.  2 


Evolution  of  self-consistent  equations  for  spectrum 
and  field  perturbations. 


Electron  distribution  and  spectrum 


Electron  and  ion  density  profiles 


Equilibrium  and  quasi-equilibrium  electric  fields 


Quasi-equilibrium  condition 


Fig.  3  -  The  Gaussian-Maxwellian  magnetoplasma  column. 
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13  ABSTRACT 


The  problem  of  wave  dispersion  and  stability  for  a  class  of  hot,  inhomogeneous, 
collisionless  magnetoplasmas  is  reduced  to  the  solution  of  an  integral  equation 
with  well-behaved  kernel.  Admissible  configurations  include  those  for  which  the 
externally  applied  and  internal  ambipolar  fields  form  a  generalized  harmonic 
oscillator.  The  full  s,  t  of  Maxwell's  equations  is  used  to  arrive  at  self- 
consistent  perturbation  fields  in  terms  of  the  equilibrium  particle  distributions. 

An  illustrative  example  treats  a  magnetoplasma  column  with  Gaussian  radial  profile 
and  Maxwellian  velocity  distribution  in  a  state  of  quasi-equilibrium. 
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